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ABSTRACT 

The primordial non-Gaussianity signal, if measured accurately, will allow us to distinguish between 
different candidate models for the cosmic inflation. Since the galaxy groups located in void regions 
are rare events, their abundance may be a sensitive probe of the primordial non-Gaussianity. We 
construct an analytic model for the mass function of void groups in the framework of the extended 
Press-Schechter theory with non-Gaussian initial conditions and investigate how it depends on the 
primordial non-Gaussianity parameter. A feasibility study is conducted by fitting the analytic mass 
function of void groups to the observational results from the galaxy group catalog of the Sloan Digital 
Sky Survey Data Release 4 with adjusting the primordial non-Gaussianity parameter. 
Subject headings: cosmology: theory — large-scale structure of universe 



1. INTRODUCTION 

In precision cosmology, a fundamental task is to distin- 
guish between various candidate models for the cosmic 
inflation. One way to perform this task is to measure 
the degree of non-Gaussianity in the primordial den- 
sity field. Although the primordial density field is re- 
arded as nearl y Gaussian in all inflationary scenarios 
Guth fc Pi|[l982 L the de gree of its non-Gaussianity dif- 
fers between the models. For instance, in the single-field 
slow-roll inflationary model, the deviation from Gaus- 
sianity is small enough to be unobservable (jVerde et al.l 
120011) . Meanwhile in some multi-field inflationary mod- 
els the non-Gaussianity can be generated to a detectable 
level (Battefeld & Easthcr 2007, and references therein). 
For a thorough review on the predictions of various in- 
flationary models for the primordial non-Gaussianity, see 
iBartolo et aP (12004 ) 

The cosmic microwave background radiation (CMB) 
is a useful probe of the primordial non-Gaussianity 
as it reflects the linear density perturbations. Yet, 
if the primordial non-Gaussianity is scale-dependent, 
the constraint from the CMB an alysis is rather re- 
stricted to the ve ry large-scale ()Verde et al.l 120011 : 

iKomatsu fc Spergeil I2001D . The large-scale structure 
is a powerful alternative probe of the primordial 
non-Gaussianity on the sub-CMB scale. It has been 
well known that the abundance of high-z clusters 
are rare eno ugh to constrain th e primordial non- 
Gaussianitv (iLucchin fc Matarres^ 119881: IChiu et al 



29981: iMatarrese et al. 2000': We inberg fc Kamiqnkowski 



20031: ICarbone et al.. 
Grossi et al.l I2OO90 



200& ILoVerde et all I2008I .- 
This probe, however, is likely 
to suffer from large systematics involved in the inac- 
curate measurement of the masses of high-z clusters. 
The abundance of cosmic voids is another probe of the 
primordial non-Gaussianity based o n the large-scale 
structure (jKamionkowski et al.l l2009( l. One difficulty 
in using this probe lies in the fact that ther e is n o 
unique way to define voids (jCarbone et al.l I2008D . 
Recently, it has been claimed that using the clustering 
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properties of highly biased large-scale structures the 
primordial non-Gaussianity parameter can be measured 
with accuracy as high as the one obtained from the 
CMB analysis (Carbo nTeFaLl [20081 : ISlosar et all [20081 : 
iJeong fc Ko matsu 2009). 

Here, we propose the mass function of present galaxy 
groups embedded in void regions as a new probe of the 
primordial non-Gaussianity. It is not only the high-z 
clusters but also the low-z void groups that are so rare 
that their abundance may depend sensitively on the ini- 
tial conditions. Furthermore, the mass estimation of low- 
z galaxy groups sh ould be much mor e reliable than that 
of high-z clusters (jYang et al.ll2007D . T hroughout this 
Lette r, we assume a WMAP 5 cosmology (jDunkley et al.l 
[20091) . 

2. AN ANALYTIC MODEL 

The Press-Schechter theory (i Press fc Schechted 11974 
PS hereafter) provides an analytic framework within 
which the number density of bound objects as a func- 
tion of mass, dNps/dM, can be obtained: 



dN _ p d 
dM ^ MdM 



p{dM)dS 



(1) 



where p is the mean background density, 5c{z) is the 
critical density contrast for gravitational collapse at red- 
shift z, p{Sm) is the probability density distribution of 
the density field Sm smoothed on the mass scale of M, 
and A is the normalization constant. Basically, equa- 
tion (|T]) states that the number density of bound ob- 
jects can be inferred from the differential volume frac- 
tion occupied by those regions in the linear density 
field whose average density contrast Sm reaches a cer- 
tain threshold, Sc- In the original PS theory, the ini- 
tial density field is assumed to be Gaussian as p{6m) — 
exp [-Slj/{2al,)] H^/^^gm)- 

Here gm is the rms density fluctuation smoothed on 
the mass scale M and (5c = (^coZ-D+l-z) where (5co is the 
critical density contrast at z = and -C>+(z) is the lin- 
ear growth factor. For a WMAP 5 cosmology, we find 
(5co ~ 1.62. Note that the mass function of bound objects 
depends on the initial conditions through its dependence 
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on aM which is a function of the density parameter, flm 
and the amphtude of the hnear power spectrum, (Tg. 

Now, let us consider the case of non-Gaussian initial 
conditions that is often characterized by the primor- 
dial non-Gaussianity parameter /nl as V'(x) — 0(x) -|- 
/jVL ((^^(x) — (0^(x))) where is a Gaussian random 
field and V'lx) is the li nearly extrapolated gravitational 
poten tial at z = (e.g.. lLoVerde et 311120081 iGrossi et al.1 
I2009D . The functional form of p{5m) for the non- 
Gaussian case will be in general different from the Gaus- 
sian case. However, provided that the degree of the non- 
Gaussianity is very small and scale independent, p{5m) 
has the same f unctional form as the Gaussian case at first 
order ( Lucchi n fc MatarresellT988l: iMatarrese et al.ll2000l : 
IVerde et al.li200ll) . The onlv difference is the value of the 
critical density contrast, 6c* , which is related to that of 
the Gaussian case, Sr as 



S„{z) = Sc{z) 



Sc{z) 



1/2 



(2) 



Here S3 is a skewness parameter, related to the 
primordial non-Gaussianity parameter /nl as S3 — 

|/nlM3^V(A'2^^)^ where n'^^ and ^3^^ denote the vari- 
ance and skewness of the smoothed non-Gaussian den- 
sity field at first ord er, respectively (see eqs. 43-45 in 
IMatarrese etalll2000[ ). Therefore, the PS mass func- 
tion with non-Gaussian initial conditions is a function 
of M and /nl: dNpsihh, M)/dM. The case of /nl = 
corresponds to the original one with Gaussian initial con- 
ditions. 

On the group scale, dNpsifyh, M)/dM does not 
change sensitively with the initial conditions since the 
galaxy groups are not rare events. However, those groups 
located in voids should be so rare that their abundance 
may depend sensitively on the initial conditions. Before 
deriving the abundance of void gro ups, we cla r ify th e 
meaning of a void region. Following iHahn et al.1 (|2007f ). 
we define a void region on mass scale M as a region where 
the three eigenvalues of the tidal tensor at a given region 
on mass scale M are less than zero. Then, we replace 
p{6m) in eq. by p{Sm\X'ij^.j, < 0) which represents 
the conditional probability density distribution that the 
density contrast has a certain value on the mass scale M 
provided that the largest eigenvalue of the tidal tensor is 
negative on some larger mass scale M' > M. From here 
on, 6 denotes the density contrast on the mass scale M, 
while X'l, A2, A3 are the three eigenvalues of the tidal 
field on some larger mass scale M'. 

The joint probability distribution p(6, X[, X'2, X'3) for 
the case o f Gaussian initial conditions has been already 
derived by[Si (|2006l ) as 



piS,X[,X'2,X'3)^ 



1 



3375 



exp 
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A'2)(Ai-A^)(A'i-A^,( 



with (T^ 



<7^~ cr'2, I[ ^ X[ + A'2 + A(,, and I[ = X'^X'^ + 
A2A3 -I- X'iX3. Here a and a' represents the rms density 
fluctuations on the mass scale M and M', respectively. 
The conditional probability density, p{S\Xi < 0), is now 
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Fig. 1. — Mass function of void groups (solid line) and of all 
groups (dashed line) in the extended Press-Schechter formalism. 



written as 



p{d\X[ < 0) 



/°oo ^Ai J^i dX, J^i dX3p{S, X[, X'„ X'3) 

f_^dx,j^l^dx,f^idX3p{x[,x'„x',) ■ 

(4) 

Putting p{5\X[ < 0) in eq. ([1]), we evaluate the mass 
function of void groups with Gaussian initial conditions, 
dN^/dM. 

Usin g the methodology suggested by iLoVerde et al.l 
(l200l . we model the mass function of void groups with 
non-Gaussian initial conditions, dNy'^ /dM, as 



dN^"" _ dN^ dNpsihL,M)/dM 



dM 



dM diVps(/NL = 0,M)/dM' 



(5) 



where diVps(/NL = Q,M)/dM and dNpsihh, M)/dM 
represents the PS mass function with Gaussian and non- 
Gaussian initial conditions, respectively (see eq. [4.20] in 
Loverde et al. 2009). Figure [1] plots the mass func- 
tion of void groups (solid line), which normalized as 
/ dNv / dXogMi — Nyg where N^g is the total number 
of void groups found in the observational data (see §3.) 
The mass function of all groups (dashed line) are also 
plotted for comparison. As can be seen, the mass func- 
tion of void groups decreases very rapidly with mass, 
which indicates that it must depend sensitively on the 
initial conditions. 

It has to be mentioned here that equation ^ has 
yet to be vali dated against numer ical results. The 
methodology of ILoVerde et al.l (j2008D that equation (O 
is based on has recently been tested against N-body sim- 
ulatio ns and found to b e valid in the high-mass sec- 
tions (jGrossi et al.ll2009l ). Yet, to fully justify the use 
of equation ([5]) for the evaluation of the abundance 
of void groups with non-Gaussian initial conditions, it 
will be required to test equation ([5]) numerically on the 
group-mass scale. Furthermore, for the case of scale- 
dependent non-Gaussianity the mass function of void 
group s would have rnuch m ore complicated even at first 
order (jLoVerde et al.ll2008[ ). The focus of this work, how- 
ever, is on the proof of a concept that the abundance of 
present galaxy groups in voids can be used as a probe of 
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Fig. 2. — Mass function of void groups with the best-fit value 
of /nl (solid line) determined by fitting to the observational data 
from the SDSS DR4 (dots). The errors include both the Poisson 
noise and the cosmic variance. The mass function of void groups 
for the Gaussian case (/nl = 0) is also plotted for comparison 
(dashed). 

the primordial non-Gaussianity. Henceforth, here we use 
equation ^ as an ansatz and consider only the scale- 
independent case for simplicity. 

3. A FEASIBILITY STUDY 

We conduct a feasibility study by comparing the an- 
alytic mass function of void groups obtained in §2 with 
the observational result from the galaxy group catalog of 
the Sloan Digital Sky Survey Data Release 4 provided by 
(|Yang et al.ll2007D who measured the masses of the SDSS 
galaxy groups were estimated from the characteristic lu- 
mino sity (or ste llar masses) using the WMAP 3 cosmol- 
ogy (jSpergel et al. 2007). To identify void groups from 
the SDSS group catalog, we use the rea l-space tidal field 
reconstructed bv lLee fc Erdogdiil (|2007l ) on 64"^ pixels in a 
box of linear length 40 0/t~^Mpc from the T wo Mass Red- 
shift survey (2MRS) (" Huchra et aI1l2005D . We smooth 
the 2MRS tidal field with a Gaussian filter of scale radius 
8ft-~^Mpc on which scale the density field is still in the 
quasi-linear regime, calculate the three eigenvalues of the 
smoothed tidal field at each pixel and mark as voids those 
pixels in which all three eigenvalues are less than zero. 
A total of 550 galaxy groups at redshifts 0.01 < z < 0.04 
in mass range of 11.7 < log A//(/i~^A/0) < 13.4 are 
found in the void pixels. Binning their mass range 
in the logarithmic scale, we measure their abundance, 
dNv/d\ogM, which is renormalized to be / dNy = 1. 
To account for the cosmic variance as well as the Pois- 
son noise in the measurement of dNv/d\ogAI, we divide 
the selected void groups into 6 subsamples and measure 
dNv/d\ogM for each subsample separately. We calcu- 
late the jackknife errors as the standard deviation in the 
measurement of the mean dNy /dlogM averaged over 
the 6 subsamples at each mass bin. 

We fit the observational result to the analytic model 
by adjusting the value of /nl- To account for the cor- 
relations between the mass bins, we employ the general- 
ized x^-statistics to determine the best- fit value of /nl: 
= [n,-n(logA//,;/NL)]Ci^^[n,-n(logAf,;/NL)] where 



Fig. 3. — Degeneracy curves in the /NL-f^m plane with the value 
of erg set at 0.76 (solid)and at 0.83 (dashed). 

Hi = dN/dlogMi and 7i(log Af^; /nl) represents the ob- 
servational and analytical results evaluated at the j-th 
logarithmic mass bin, logA^^, respectively. And (Cy) is 
the covariance matrix defined as = ((n^ — noi){nj — 
nQj)) where riQi represents the mean of rii averaged over 
all samples. Finally, the uncertainty in the measurement 
of /nl is calculated as the curvature of the function 
at the minimum. Through this fitting procedure, we find 
/nl = 36 ± 1 . Figure [2] plots the observational results 
(dots) and the analytic model with the best-fit value of 
/nl = 36 (solid line). The median redshift of the SDSS 
void groups, z — 0.03, is used for the value of z in the an- 
alytic model. For comparison, the analytic model with 
/nl = is also plotted (dotted line). As can be seen, 
the observational results agree better with the analytic 
model with non-Gaussian initial conditions. 

It is, however, expected that there is a degeneracy 
between /nl and the other key cosmological parame- 
ters on which the mass function of void groups depend. 
Here, we investigate the degeneracy between rim and 
/nl, setting erg at the value determined by WMAP 5 
cosmology. Varying the values of flm and /nl, we re- 
calculate dNv / dlogM at a typical group mass scale of 
M = IQ^^h'^MQ. Figure [3] plots a family of the degen- 
eracy curves in the r2m-/NL plane with the value of erg 
set at 0.76 (solid) and 0.83 (dashed). As can be seen, 
a strong degeneracy exists between the two parameters. 
For a fixed value of dNy / d\ogM , the value of /nl in- 
creases as the value of Vtm decreases. The comparison 
between the solid and the dashed lines indicate that the 
degree of the degeneracy between Vim and /nl increases 
as the value of erg increases. It also indicates that if f2„i 
is fixed, the value of /nl increases as the value of erg de- 
creases. To break this parameter degeneracy, it will be 
necessary to combine our analysis with other analyses. 

The preliminary results of this feasibility study, how- 
ever, are subject to several caveats. First, the analytic 
model assumes the scale-independent Gaussianity. To 
be more realistic, it is necessary to account for the scale- 
dependent non-Gaussianity. The second caveat lies in the 
limitation of the Press-Schechter formalism. As shown by 
several authors ()Lee &: Shandarinlll998t[Sheth fc TormenI 
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Il999t I Jenkins et al.|[200ll ). the real gravitational process 
deviates from the spherical dynamics on which the PS 
mass fiinction is based. It has to be tested how sig- 
nificantly the deviation of gravitational collapse process 
from the spherical dynamics affects the abundance of 
void groups. The third caveat comes from the fact that 
the val idity of e quation ^ has yet to be confirmed. Al- 
though lGrossi et al. ( 2009) have shown th at this method- 
ology suggested by iLoVerde et al.l ()2008i ) to model de- 
partures from non-Gaussianity leads to an excellent ap- 
proximation on the cluster scale, it has to be confirmed 
by N-body simulations whether the same methodology 
can be used to count the number of void groups with 
non-Gaussian initial conditions. Fourth, the different 
mass-to-light ratios of the void galaxies from that of 
the wall g alaxies has t o be taken into account. Ac- 
cording to iRoias et al.l ()2005f ) , the specific star forma- 
tion rate in void galaxies is higher than that in wall 
galaxies, suggesting that the mass of the void groups in 
the SDSS Gal axy group catalog are Hkely to be overesti- 
mated. Fifth. lYang et al.l (|2007l ) measured the masses of 
SDSS galaxy groups assuming the old WMAP 3 cosmol- 
ogy llpergcl et al. 2007). It will be necessary to use the 
values of the most updated WMAP 5 parameters for the 
more accurate calculation of the masses of void groups. 

As a final conclusion, we have proved that the abun- 
dance of void groups can in principle be a useful probe 
of the primordial non-Gaussianity parameter. For a ro- 
bust probe, however, it will be required to refine further 
the analytic model of the abundance of void groups and 
to improve the mass estimation of galaxy groups in void 
regions, which is the direction of our future work. 
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